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In this paper, we discuss the phenomenon of a spontaneous polarisation of a neutral 
hard planar interface for valence asymmetric coulombic systems. Within a field theoretical 
description, we account for the existence of non trivial charge density and electric potential 
profiles. The analysis of the phenomenon shows that the effect is related to combinatorics 
in relation with the existence of the two independent species cations and anions. This 
simple and basic feature is related to the quantum mechanical properties of the system. The 
theoretical results are compared with numerical simulations data and are shown to be in 
very good agreement, which a fortiori justifies our physical interpretation. 



I. INTRODUCTION 



There has been recently a renewed interest in the study of the structure of the double layer. This 
stems both from the existence of new important domains of application as for instance those related 
with the nanotechnologies, microfluidics, microbatteries, and electrochemical sensors, biology or 
also the use of new ionic liquid solvents [l|, 2] . These new domains represent a theoretical challenge 
as the behaviours of the charged systems does not in general follow the popular predictions of the 
Poisson-Boltzmann approach. There can be many reasons for parting from the standard behaviour. 
For instance, it can be related to excluded volume like packing effects at high local concentrations 
at the interface but also, in the regime of low reduced temperatures, to the existence of 
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strong Coulomb interactions which require to account for correlations beyond mean field theories 
like the popular Gouy-Chapman (GC) approach. In the latter case, there has been a series of 

13l | devoted to the challenging problem of the anomalous 
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behaviour of the electric capacitance. 

In this paper, we investigate another type of systems which also depart from the standard GC 
app roach. It is the case of the asymmetric in valence electrolytes. In an early paper, Torrie et al. 



14l ] have shown the existence of a polarization of the interface even at the point of zero charge 



(PZC) for an extended restricted primitive model of ions with different valences. This phenomenon 
is far from obvious and standard intuitive approximations like the GC theory are unable to describe 
it. In contrast to the case of asymmetric in size ions, where the smallest ion by coming closer to the 
interface induces a polarization of the interface, the behaviour of valence asymmetric ions is rather 
non intuitive. Torrie et al. have used a modified Poisson-Boltzmann approach to describe this effect 
which is in our opinion rather costly in its mathematical application. More recently Henderson 
et al. [l^ have addressed the problem using standard approximations in the liquid state theory 
searching for simple analytic expressions. They were able to reproduce the polarization effect for 
the size asymmetry but not for the valence asymmetry. From a practical point of view, the interest 
of studying the behaviour of the PZC is that this quantity is often considered as an indication 
of chemical interaction of the ions with the electrode so called "specific adsorption". Here, by 
considering a hard wall, we discard the role of "specific interactions" and intend to propose a 
simple physical interpretation. 

The article is organized as follows. In the first Section, we present the field theory formalism for 
valence asymmetric electrolytes and introduce the meaningful physical parameters for the system. 
In the following two Sections IIIII and llVl we derive respectively the expressions of the charge profile 
and of the electric potential accross this interface at a neutral hard wall. Then in Section |Vl we 
compare the results to the exact relation given by the charge contact theorem {isl and verify the 
consistency of our expressions. Finally in Section IVll we discuss the different profiles with respect 
to numerical simulations results [3]. 
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II. FIELD THEORY FOR POINT IONS. 



Formalism 



We consider a field theoretical description of a system of point ions situated in a half space, 
in contact with a hard wall. The dielectric constant e = £^£0 is uniform throughout the space. 
As in [ill . [1^ . we introduce a field theoretical expression of the grand potential in terms of the 
fields /3+(r) and /J-(r) representing in space the density fields for cation and anion distributions, 
respectively 



Pp±(r)exp{-/3i/[p±]} 



(1) 



where P = \/{kBT) is the inverse temperature. The grand potential is then obtained as P{—pV + 
7 A) = — ln6[/9±], where p is the pressure, V is the volume, 7 is the surface tension, and A is the 
area of the electrodes. For the study of the Coulomb interactions, it is convenient to use the charge 
density field g(r) = z+/9-|-(r) — 2;_/f3_(r), with 2+, z_ the valences of cations and anions respectively, 
and the total density field s(r) = p+{r) + P-(r). The Hamiltonian as a functional of these fields is 



PH[q, s] = s] + /?//^°"'[(?] - j I3iiss{r)dr - j p^igq{r)dr 



(2) 



where Hs = /^+ + /i-, fJ-q = -2^+^+ — and /x± are the chemical potentials of the ions. The 

entropic part of the Hamiltonian which accounts for the quantum mechanical degrees of freedom 
in the phase space jiol and the terms related to the chemical potentials can be written as 



PH'^^'ig, s\- j Pnss{v)dv - j Pfigq{r)dr = j 



q{r) + s(r) 



In 



+ 



s(r) - q{r) 



In 



q{r) + s(r) 
s(r) - q{r) 



1 



dr 



2p- 



1 



dr (3) 



where p± = exp(— /3/i±)/A^ with A being the de Broglie wavelength. The second term in the 
Hamiltonian is the coulombic contribution 

-2 qir)q{r') 



PH 



Coul 



Pe' f 



-drdr' 



(4) 



where e is the elementary electric charge. 

In order to perform the calculations we expand the Hamiltonian of equation ([2|) around the 
mean field profiles ^ = and s = p which minimize the Hamiltonian in the case of the neutral 
hard wall [ig], where p = pj^ + p-. Beyond the mean field solution, introducing 5s = {s — p)/p and 



q/p, the Hamiltonian is 

PH = pV + 



5q\v) 



dr + PSH + PH' 



Coul r 



(5) 
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where PSH contains terms of the expansion of orders higher than quadratic. We introduce Zis = 
^JzyzZ and z^s = (2+ — z_)/ yjzj^z_. The first coefficient is related to the ionic strength as 
+ z'^p_ = However this parameter is not sufficient, the same for a 2:2 or a 

4:1 electrolyte and it is the second coefficient Zas which is really characteristic of the asymmetry 
in valence between ions. With these notations we have 

(55H = ^ f [5s(v) - 5s^{r)\5q^{r)dr - ^ [ [l - 25s{r) + 56s'^{r)] 6q^{r)dr 

IS IS 

+ 4^ /(I + ^as) [1 - S-^Kr) + QSs\t)] 5q\v)dv + ... (6) 

Note that the quadratic term in the Hamiltonian remains diagonal for the q and s fields also for 
the valence asymmetric systems and odd terms in q appear only in f35H, the smallest odd power 
in q being three. 

As in [20], the Hamiltonian can be further simplified by scaling the charge density field 5q — > 
6Q = bqjzis and by defining a new unit charge e = Zigt. In this case for the even terms we recover 
the expansion of the symmetric 1:1 electrolyte and the Hamiltonian is 

PH = pV + ^J [5s\r) + 6Q\r)] dr + p6H + pH^^^'^Qiv)] (7) 



with 



P5H = -f / Hr)6QHv)dv + ^ I 5Q''{r)dr-^ J SQ\r)dr + J 6Q\r)dr + ... {8) 

where we have kept from equation ([6]), for clarity, only the terms of interest. The Coulomb 
interaction term is also modified 

Svrp J \r — r| 

where we have introduced the inverse Debye length, = (/3(2;^p+ + z^p_)e^/e)^/^ = (Ppe^ /e)^^"^ , 
which shows that has the same expression as for the 1:1 electrolyte except that the elementary 
charge is now e. Note that the scaling introduced in this section can be considered as a simple 
renormalisation of the electrostatic quantities in relation with the ionic strength. In the following, 
we perform all calculations in terms of the charge density field 5Q and transform the result in the 
physical field 6q at the end of the calculations. 

The parameter Zis has been completely absorbed in the new field 5Q and the charge e and like in 
[20I we can use the results for the 1:1 symmetric electrolyte and we mainly need to focus on the new 
terms in comparison to the symmetric case which are associated with the asymmetry parameter 
Zas- For the interfacial properties, at the one loop level of the calculations we will assimilate the 
activity p and the average density p. 
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III. CHARGE PROFILE AT THE PZC 

The calculation of the charge profile illustrates the peculiarity of the asymmetric systems. For 
the symmetric z : z electrolytes at the neutral interface, for reasons of symmetry, average quantities 
which include an odd number of charge densities, as for instance the charge profile, vanish. Another 
way of stating this is by noting that it is impossible by applying the Wick theorem to pair the 
charge fields and simultaneously calculate a quantity with an odd number of charge fields, because 
this would imply the existence of odd coupling constants. Such odd terms exist for the asymmetric 
systems as a consequence of the expansion of the entropic contribution to the Hamiltonian. The 
first contribution is due to the three body coupling in equation dH) and corresponds to the diagram 
in figure [J where we have the standard 1/2 symmetry coefficient associated and the analytic 
expression is 

(5Q(r)) = jdr' < 6Q{r)6Q{r') >< 6Q{r')6Q{r') > (10) 
where the inhomogeneous correlation functions are those given in [20|], where 

{SQir)5Qir')) = - ['^(-0"^^-^ + J e "^^^^^^^^^ 

(11) 

where distances of the points from the wall and R, R' are respectively the projections 

of r and r' parallel to the wall and where K' = ^ K"^ + K'j-^ . In the previous correlation, when 
calculated at the same point, we renormalize the dirac distribution as indicated in [l^. The bulk 
part which does not depend on the distance to the wall, vanishes with the rest of the graph as 
a consequence of the electroneutrality condition J {6Q{y)6Q{y')) dr' = 0. And finally the excess 
contribution at the wall gives jlj] 

{SQ{r)SQ{r)) = ^_I{x) (12) 



where we used the reduced distance to the wall x = xKd and I{x) = e "^^^/{t + — l)'^dt. 
We finally obtain our result, the profile for the charge 

q{x) = ^ZasZis rjp F{x) (13) 

where the coefficient rj = K'^/{87rp) and the function F(x) = —21{x) + f{x) + cq where 

dt 



1 [i2t-l) (2i-l)(2i + l) 



(t + x/t^)2 



(14) 
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and Co = -7r^/3/12 + In 3/4 + 1/4 w 0.0712. 

The generic form of the function F(x) is given in figure [21 At the wall, we have -F(O) = 
— 2/3 + 2co, thus the charge density at contact is 

q{0) = ZasZisPV ^co - ^ (15) 

We also note that the integral of F is zero implying that as expected the profile verifies the 
electroneutrality condition J q{x)dx = 0. 

IV. POTENTIAL PROFILE AT THE PZC 

The electric potential can be obtained using the charge profile from the exact expression 

(5eil]{x) = -— / {x' - x)q{x')dx' (16) 

^ J X 

using equation (fT3]) . we obtain 

^eV'(x) = - — ^ r {x' - x)F{x')dx' (17) 
= '-^vG{x) (18) 



where 



G{x) = / ' . cie-^ (19) 



with ci = 7r^/3/24 + In 3/8 - 1/4 0.11405. The generic form of the function G{x) is given in 
figure [3l We note that the derivative of this electric potential, which is proportional to the electric 
field, is zero at the interface as expected for a neutral hard wall. And at the contact with the wall, 
we have 

/3eV'(0) = -— ??co (20) 

Zis 

V. CHARGE CONTACT THEOREM 

In [3], we have derived a charge contact theorem which relates the contact value of the charge 
profile to the electric field accross the interface. Using this exact relation, in the following, we verify 
the consistency of our expressions for the charge density and for the electric potential profiles. For 
an asymmetric in valence system of point ions, the charge contact theorem reads 

g(0) = Pej^{zlp+{x) + zlp^{x)) (^|^) dx + f3{z+P+ + z_P_). (21) 



Note that in comparison to the notations in [l^, here, g is a density and does not include the 
electric charge. At the lowest order in the loop expansion, we can take the profile of the density 
and the electric potential respectively at the zero loop order which corresponds to a constant profile 
and at the first order, which is our result equation ([TH]) . So the first term becomes 



- Pezl,pip{0) = (z+ - z^)prico = ZasZispijco (22) 
For the second term, at the same order, the partial pressures are given by the Debye approximation 

^P. = ^ (23) 

z+ + Z- 

where l3P = -K%/{2A-k) = -pr]/^ is the total pressure. We thus have 

(3{z+P+ + z.P.) = (3{z+ - z^)P (24) 
= -^ZasZis-qp (25) 

We can see that the sum of the two contributions equations (|22p and (|25p gives the value obtained 
from the direct calculation of the charge density profile given in equation (jl5|) . This shows the 
consistency of our results for the charge density and electric potential profiles obtained at the first 
order in the loop expansion with the exact contact charge relation. 



VI. DISCUSSION 

In section IIIII and IIVI we have shown the existence of a spontaneous polarisation of a neutral 
interface. This phenomenon is absent in symmetric systems where we have the same number of 
anions and cations at each point accross the neutral interface, although we have shown there is a 
profile for the total number of ions [3]. The polarisation is then directly related to the asymmetry 
in valence of the ions, that is to the deplacement of the equilibrium in number of the ionic species 
in order to satisfy the electroneutrality condition. A way of understanding this profile is to consider 
the depletion of the ionic profiles due to the electrostatics existing at interfaces [IS] . This depletion 
is stronger and shorter ranged the higher the charge of the ions. Assuming z+ > for the more 
highly charged cations, we understand that their profile is more depleted closer to the interface 
leading to a negative charge profile at contact. This is followed by a distribution of negative ions 
which is then electrically compensated by positive ions further away from the wall. The electric 
potential at the wall is then negative and increases to reach zero towards the bulk. Such qualitative 
statements have been discussed on the occasion of the contact theorem for asymmetric electrolytes 
211. 
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TABLE I: Extrapolation of the PZC by LLSQ fit for the 2:1 electrolyte from numerical data results in \L 



concentration 


PZC 


V 


0.05 Mol/1 


-0.0351 


0.909 


0.5 Mol/1 


-0.0989 


2.87 


1.0 Mol/1 


-0.130 


4.07 



More quantitatively, we have seen in Section |Vl that the profiles are such that they verify 
consistently the charge contact theorem at given order of approximation. From equation (jl3p 
and (jlSp . the expressions of the charge density and the electric potential profiles are analytic 
and simple. Here are some of their characteristic features. The two profiles depend respectively 
on the universal functions F and G in terms of the reduced distance x and linearly on the loop 
expansion parameter rj. As a consequence, we can see that similar to the desorption phenomenon 
and to the anomalous capacitance behaviour [2^ this polarisation of the interface increases as 
the reduced temperature decreases. We recall that the reduced temperature is proportional to 
the temperature and to the dielectric constant and inversely proportional to the product of the 
ionic charges. Therefore, there is a fair number of physical systems which correspond to a low 
reduced temperature as it is possible to decrease the systems' temperature or dielectric constant 
or conversely consider highly charged ionic species. Typically for a 0.05 molar solution of a 2:1 
electrolyte at room temperature in water we have rj ~ 0.909. Then, the value of the profiles are: 
linear on the difference on the product ZasZis = Z-^. — z- for the charge profile where as the reduced 
electric potential scales on the ratio Zas/zis = {z+ — z^) / {zj^z^). As a consequence the behaviour 
of the charge and of the electric potential are different. The charge scales proportionally to the 
difference 2;+ — z-^ where as the ratio Zas/ Zis is bound by 1 for high charge asymmetries. 



We have compare d q uantitatively our point ion expression for the PZC with the numerical 
simulation results in [l5l]. Note that for this comparison, the contact value distance corresponds 
to half the ionic diameter for the finite sized ions of the MC simulations, where as for the point 
ions it corresponds to a vanishing distance to the wall. As the size asymmetry is immaterial for 
the point ions, we first extrapolate the simulation results for the equal size ions using a linear least 
squares quadratic (LLSQ) fit for the 2:1 and 3:1 electrolytes as shown in figure H] and O The fits 
give for the various concentrations used in [l^ the PZCs' shown in Table I and II, were we have 
also given the value of the parameter r/. For all these systems we note that > 1, which indicates 
the importance of beyond mean field correlations. In contrast to the numerical simulation results. 
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TABLE II: Extrapolation of the PZC by LLSQ fit for the 3:1 electrolyte from numerical data results in \L 



concentration 


PZC 


V 


0.033 Mol/1 


-0.0726 


1.57 


0.333 Mol/1 


-0.226 


4.98 


0.666 Mol/1 


-0.320 


7.04 



simple, analytic approaches of the liquid state theory like the LMGC (Linear Modified Gouy 
Chapman) and MSA (Mean Spherical Approximation) theory give a vanishing value of the PZC 

The comparison of the extrapolated values of the potential with the field theory are given as a 
function of the density in figure [6] and [71 For nearly all values, the results of the field theory for 
the simple point ions are very close to those of the simulations, for the concentrations investigated 
and for both 2:1 and 3:1 electrolyte. The only point which is not as accurate corresponds to the 
highest density considered for the 2:1 electrolyte. The overall agreement of the theory seems to 
indicate, that the excluded volume effects are not crucial. And favours of our interpretation that 
this polarisation phenomenon is essentially entropic in nature meaning it is related to the number 
balance of the ionic species, which follows the valence asymmetry. 



Given the results of the simple point ion model, we have considered extending the comparison 



15(. In this article, 



for non symmetric electrolyte both in size and in valence as they are studied in 
for valence symmetric systems the LMGC and MSA approaches give a reasonable dependence 
with respect to the ionic diameter ratio. It is then tempting to assume an additive approximation 
combining the size asymmetry well described by standard liquid state theory approaches and 
the valence asymmetry which is well described by our field theory approach. We have used this 
assumption for the results presented in figures [8] - [TTJ 

In figures [5] and El the theory is obtained by adding the PZC values of the LMGC approach to 
those of the field theory (Tables I and II), respectively for the 2:1 and 3:1 electrolytes and for the 
corresponding concentrations. For both types of electrolytes, this approximation gives good results 
for near symmetric in size systems. At the intermediate density considered, the agreement extends 
over all diameter ratios whereas it respectively underestimates and overestimates the effect at low 
and high densities. 

In figures [TO] and [TTl we replace the PZC values of the LMGC by those of the MSA approach. In 
this case, the results are overall much closer to the simulation results for both type of valences and 
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all concentrations. However, in the case of the 2:1 electrolyte the approximation is more succesful 
at lower concentration. Whereas for the 3:1 electrolyte, the comparison is incredibly good over the 
whole concentration range considered. 

Thus simply adding the PZC effects due to the size and valence asymmetries gives rather good 
results and allows with a simple model the simultaneous description of both asymmetries. We 
believe that the discrepancy appears to be lower for the 3:1 electrolyte as in this case the valence 
effect which is rather well predicted by the field theory is larger and possibly may partly hide the 
discrepancies in the description of the size effects. 
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VII. CONCLUSION 

In this paper, we have studied the behaviour of valence asymmetric electrolytes at simple hard 
neutral interfaces. We have discussed the existence of a polarization of the neutral interface for 
these systems which is not related to any specific interaction with the interface. 

Using a point ion model, we show that this effect is purely entropic. More precisely, it is not 
related to volume exclusion potential effects rather it is connected to the number balance between 
the ions. This phenomenon is quite fundamental as its origin comes from the combinatorics and 
counting of the states in the phase space set by basic rules of the quantum mechanics. These familiar 
rules give the entropic contribution of the particles of the same kind which are indiscernible and as 
a consequence also set the combinatory for anions and cations which are discernible species. The 
physics of these rules appears in the entropic part of the functional of our field theory Hamiltonian. 
The role of this entropy has previously been emphasized discussing the desorption phenomenon 
and the related anomalous capacitance behaviour. In this case, the entropic functional leads to a 
coupling between charge and total density fields 
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22l ]. In the present work, the same entropic 



functional gives couplings for the charge field alone and gives the correction to the trivial vanishing 
mean field charge density profile. Another interesting aspect, in the approach, is that on the field 
variables we are allowed algebraic operations (e.g. linear combinations of the fields like the charge 
or the total density) which turn out to be meaningful physically and are less straightforward to 
implement when discussing particles. Finally, for the charge and the electric potential profile, we 
obtain expressions which are analytic and rather simple with well identified scaling behaviours in 
terms of the physical parameters. In particular, we discuss the role of two parameters related to 
multivalency. One is the product of the valences, which is related to the ionic strength in the 
system. The second is more characteristic of the asymmetry as it is proportional to the difference 
of the valences. 

We believe that this study illustrates the interest of our FT approach. One important aspect 
is the original description of the quantum mechanical degrees of freedom which leads to new inter- 
pretations and allows for certain systems to obtain a clear physical picture and simple application 
in comparison to other approaches. 
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FIG. 1: Diagram for the calculation of the charge density profile as described in the text. 




FIG. 2: F function, proportional to the charge density profile, as a function of the distance to the wall in 
reduced units. 




FIG. 3: G function, proportional to the electric potential profile, as a function of the distance to the wall in 
reduced units. 
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FIG. 4: PZC data from [l5| (circles, diamonds and squares) for the size asymmetric 2:1 electrolytes for 
concentrations 0.05, 0.5, 1.0 Mol/1 respectively extrapolated to the equal diameter systems (empty symbols). 
The dotted lines are the LLSQ fit used. 
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FIG. 5: Same caption as in figure |4] for the 3:1 electrolyte and concentrations 0.033, 0.333, 1.0 Mol/1. 
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FIG. 6: PZC for the extrapolated MC simulation results (empty circles) of figure H] and the field theory 
Table I (crosses) for a 2:1 electrolyte as a function of the ionic concentration. 
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FIG. 7: PZC for the extrapolated MC simulation results (empty circles) of figure [5] and the field theory 
Table II (crosses) for a 3:1 electrolyte as a function of the ionic concentration. 
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FIG. 8: PZC as a function of the diameter ratio, for a 2:1 electrolyte and the same physical parameters 
temperature, dielectric constant as in [15]. The symbols (circles, diamonds and squares) are for the numerical 
simulations data and the empty symbol for the extrapolated value for equal diameters, and the curves (full, 
dotted-dashcd and dashed) are for the field theory plus the LMGC theory [l^ presented in the order of the 
concentrations 0.05, 0.5, 1.0 Mol/1. 



!1 




FIG. 9: PZC as a function of the diameter ratio, for a 3:1 electrolyte and the same physical parameters 
temperature, dielectric constant as in [l^. The meaning of the symbols and curves are those of figure [8] 
except that the concentrations are 0.033, 0.333, 1.0 Mol/1. 




FIG. 10: PZC as a function of the diameter ratio, for a 2:1 electrolyte, with caption identical to that in 
figure [8] except that the MSA approximation is used in place of the LMGC approximation. 
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FIG. 11: PZC as a function of the diameter ratio, for a 3:1 electrolyte, with caption identical to that in 
figure [9] except that the MSA approximation is used in place of the LMGC approximation. 



